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ABSTRACT 

Using a method to discover and classify supernovae (SNe) in galaxy spectra, we find 
90 Type la SNe (SNe la) and 10 Type IIP SNe among the - 700,000 galaxy spectra in 
the Sloan Digital Sky Survey Data Release 7 that have vESPA-derived star-formation 
histories (SFHs). We use the SN la sample to measure SN la rates per unit stellar mass. 
We confirm, at the median redshift of the sample, z = 0.1, the in verse dependence on 
galaxy mass of the SN la rate per unit mass, previously reported by Li et al. (2011b ) for 
a local sample. We further confirm, following Kistler et al.| ( 2011 1, that this relation can 
be explained by the combination of galaxy 'downsizing' and a power-law delay-time 
distribution (DTD; the distribution of times that elapse between a hypothetical burst 
of star formation and the subsequent SN la explosions) with an index of —1, inherent to 



the double-degenerate progenitor scenario. We use the method of Maoz et al. (20111 



to recover the DTD by comparing the number of SNe la hosted by each galaxy m 
our sample with the vESPA-derived SFH of the stellar population within the spectral 
aperture. In this galaxy sample, which is dominated by old and massive galaxies, we 
recover a 'delayed' component to the DTD of 4.5 ± 0.6 x 10" 14 SNe M Q _1 yr" 1 for 
delays in the range > 2.4 Gyr. The mass-normalized SN la rate, averaged over all 
masses and redshifts in our galaxy sample, is Ris.,m{z = 0.1) = 0.10±0.01 SNuM, and 
the volumetric rate is Ruy{z = 0.1) = 0.25 ± 0.03 x 10~ 4 SNe yr" 1 Mpc~ 3 . This is 
the most precise SN la rate measurement at this redshift, and is consistent with rates 
and the rate evolution from other SN la surveys. 
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1 INTRODUCTION 

The nature of the stellar system that ends up exploding as a 
type la supernova (SN la) is still uncertain (see |Howell|20 1 1| 
and |Maoz fc Mannucci|[2011| for recent reviews). While the 
progenitor is most probably a carbon-oxygen white dwarf 
(CO WD; |Nugent et al.|[20lT] |Bloom et al.||2012[ ), in order 
for it to explode it must be ignited. The current consen- 
sus is that the WD is in a binary system, through which 
it accretes mass until the pressure or the temperature in 
the WD core rise high enough to ignite carbon and initi- 
ate a thermonuclear runaway. The two leading scenarios for 
the nature of the progenitor system are the single degener- 
ate scenario (SD; |Whelan fe Iben||1973| in which the WD 
accretes mass from a main-sequence, helium, or giant star, 
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while in the double degenerate scenario (DD; Iben & Tu- 
tukov||T984"l |Webbink||T984"l ) the WD merges with a second 
CO WD through loss of energy and angular momentum to 
gravitational waves. 

Theoretical and observational evidence exists both for, 
and against, each of the progenitor scenarios. The SD sce- 
nario, long held to be the more likely progenitor option, 
has been called into question by some recent observations. 
Pre-explosion images and early multi-wavelength data for 
SN2011fe have ruled out a red giant as the binary compan- 



ion in that explosion (Li et al. 2011c Chomiuk et al. 2012 



Horesh et al.|2012| [Brown et al.|2012f . All possible SD com 
panions, including main sequence stars, have been ruled out 
in the case of the SN la remnant SNR 0509-67.5 in the Large 



Magellanic Cloud (Schaefer & Pagnotta 2012 but see Di Ste 



fano fc Kilic|2012|and|Shappee, Kochanek, fc Stanek|20T2 

Sternberg et al. (20111 have made a statistical argument 
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that up to 25 per cent of SNe la in spiral galaxies may have 



sodium absorption lines in their spectra (see also Patat et al 



2007 Simon et al. 2009b, the result of interaction with cir- 



cumstellar material which could originate from a wind from 



a SD companion. |Dilday et aL (20121 have found evidence 
of circumstellar material around the SN la PTFllkx, which 
they explain as the result of a symbiotic nova progenitor, 
similar to RS Ophiuchi. However, such events probably con- 
stitute only 0.1-1 per cent of all SNe la. 

As for the DD scenario, |Badenes fc Mao"z| ( |2012| ) have 
recently measured the Galactic merger rate of WD-WD bi- 



nary systems in the Sloan Digital Sky Survey (SDSS; York 
|et al.|2000| | and found it to be similar to the Galactic SN la 
rate (although only ~ 10 per cent of the merger rate involves 
mergers above the Chandrasekhar mass). The DD scenario 
has long been disfavored by theoretical studies that predict 
that the merger of WDs of unequal mass would lead to an 
accretion-induced collapse to a neutron star, rather than to 
a SN la (jNomoto fc Iben][T985l |Shen et aL|[2012] see, how- 



ever, |Pakmor et al.|2010||2011[|2012| | van Kerkwijk, Chang, | 



fc Justham|2010| >. 

Some evidence favoring the DD scenario has emerged 
from reconstructions of the delay-time distribution (DTD; 
the distribution of times that elapse between a hypothetical 
5-function-like burst of star formation and the subsequent 
SN la explosions). Different progenitor models predict dif- 
ferent forms for the DTD (see |Wang fc Han|2012| for a recent 
review). A power-law DTD with an index of ~ —1 is inherent 
to the DD scenario, if the distribution of initial separations 
between the merging WDs is roughly a power law of that 
index (see, e.g., |Maoz fc Mannucci||2011[ ). DTDs from SD 
models generally cut off after a few Gyr, although extended 
DTDs in the SD context have been predicted as well (e.g., 



Di Stefano, Voss, fc Clacys 2011; Hachisu, Kato, fc Nomoto| 



2012) 



Recent studies, using various SN samples, enviro 
ments, and redshift ranges, have c onsistently measured a 



power-law DTD of the form t~ ^ (see Maoz & Mannucci 2011 
for a review). A measurement of the SN la rates in ellipti- 
cal galaxies as a function of luminosity- weighted galaxy age 



by Totani et al. ( 2008 \ was best-fitted by a power-law DTD 
of the approximate t~ x form. Galaxy clusters, where most 
of the star formation is thought to have occurred in one 
burst at z ~ 3, have also yielded SN la rates consistent with 
such a power-law DTD ( Sharon et al.||2010 Maoz, Sharon, 



|fc Gal-Yam||2010| [Barbary et al.||2012b| |Sand et al.||2012 

Measurements of the volumetric rate of SNe la in field galax- 
ies as a function of the c osmic star formati on history (SFH) 
point to a i _1 DTD [see Graur et al.|2011 for a compilation 



of such surveys and their results up to 2011, and |Perrett| 
et al.|2012| (P12) for subsequent surveys out to z ~ 1]. 

The dependence of SN la rates on host-galaxy mass and 
star-formation rate provides a further probe of the DTD. 



Mannucci et al. ( 2005 1 and Mannucci, Delia Valle, & Pana- 
gia| ( |2006[ ) showed that a broad distribution of delay times is 
required in order to explain the increase in the SN la rate per 
unit stellar mass in bluer galaxies, with at -1 DTD providing 
an excellent fit (see |Maoz fc Mannucci|2011| ). |Li et al.| ( |2011b| 
Lll) measured SN la rates per unit stellar mass in the local 
Universe and found that they followed a 'rate-size' relation, 
whereby less-massive galaxies have higher mass-normalized 



showed that the rate-size relation could be explained by the 
combination of galaxy 'downsizing' (i.e., older galaxies tend 
to be more massive than younger galaxies) and a power-law 
DTD of index -1. 



Brandt et al. (20101 and Maoz et al. (2011 Mil) intro- 



duced a method for recovering the DTD while accounting 
for the detailed SFHs of the individual monitored galaxies. 
Mil used this method on a subsample of the Lick Observa- 
tory SN Search sample (LOSS; [Leaman et al.|20lT Li et al. 
|2011a| Lll), from which they recovered a DTD with both 
'prompt' (t < 420 Myr) and 'delayed' (t > 2.4 Gyr) com- 
ponents. In |Maoz, Mannucci, fc Brandt| ( |2012| M1 2), this 
method was applied to the SDSS-II sample of SNe ( |Dilday| 
et al.|2010||Sako et al.|2011[ ), from which a similar DTD with 
prompt, delayed, and 'intermediate' (0.42 < t < 2.4 Gyr) 
components was recovered, consistent with the i _1 form of 
previous studies. 

Most SN surveys have been based on imaging, whether 
by monitoring specific galaxy samples, or volumes of space 
and the galaxies that they include. However, it is also pos- 
sible to discover SNe serendipituously among large samples 
of galaxy spectra, in which the galaxy region covered by 
the spectral aperture happens to host a SN during the time 
of exposure. Advantages of such samples are that the spec- 
tra, and hence classification, of transients is immediately 
available. Furthermore, the availability of the spectra of the 
stellar populations in the specific galactic regions that were 
effectively 'monitored' permits deriving accurate rates and 
DTDs. 

In this paper, we discover a new SN sample by search- 
ing for, and classifying, transients (including, but not lim- 
ited to, SNe la) among ~ 700,000 galaxy spectra from the 



7th data release of the SDSS (SDSS DR7; Abazajian et al 



2009[ |. |Madgwick et al.] ( j2003| M03) were the first to discover 
SNe la in SDSS galaxy spectra (in SDSS DR1). They found 
19 SNe la among ~ 100,000 galaxy spectra, from which they 
estimated the SN la rate at z ~ 0.1. Krughoff et al. (2011 
Kll) made a similar measurement, using the larger galaxy 
sample in SDSS DR5, finding 52 SNe la among ~ 350,000 
galaxies. As detailed below, in this paper we analyse a larger 
sample, using a number of improvements in transient detec- 
tion, classification, and rate analysis. We provide the details 
of our galaxy sample in Section [2] In Section [3] we describe 
our discovery and classification routine and assess its effi- 
ciency and purity. In Section [4] we present our SN sample. 
Using our SN la sample, in Section [5] we measure mass- 
normalized SN la rates, convert them to a volumetric rate, 
and compare them to previous results. In this section we 
also recover the delayed component of the DTD using the 
direct recovery method of Mil. We summarize and discuss 
our results in Section [6] 

Throughout this paper we assume a A-cold-dark-matter 
cosmological model with parameters J1a = 0.7, fi m = 0.3, 
and Ho = 70 km s~ Mpc~ . Magnitudes are on the AB 
system ( Oke fc Gunn|1983 l, unless noted otherwise. 



SN la rates than more massive galaxies. Kistler et al. (2011 1 



2 GALAXY SAMPLE 

For the purpose of this work we used the 776,447 galaxy 
spectra in SDSS DR7 which have SFHs obtained by |Tojeiro] 
et al.M2009| T09) using the VErsatile SPectral Analysis code 
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(vesp^ [Tojeiro et alj]2007| >. As in Mil and M12, we use 
the total formed stellar masses fit by vespa to the SDSS 
spectra in three time bins: 0-0.42 Gyr, 0.42-2.4 Gyr, and 
> 2.4 Gyr (corresponding to bins 24 + 25, 26, and 27 in 
T09). We specifically use the vespa reconstruction that as- 
sumes a single dust component and utilises t he |Maraston| 
( 2005 1 spectral synthesis models. A Chabrier ( 2003 1 stellar 



initial mass function is assumed in T09, and this assump- 
tion therefore propagates into the mass normalization of the 
SN la rates we report. T09 measured the total mass in the 
SDSS fiber aperture and then applied an aperture correc- 
tion to obtain the total mass in the entire galaxy (see their 
equation 22). Since, for rate purposes, we are interested in 
the stellar mass inside the fiber aperture (which was the 
stellar mass that was effectively monitored for SNe), we re- 
verse this aperture correction, and propagate the uncertain- 
ties accordingly. We reduce the masses in each time bin by 
a further factor of 0.55, due to the different flux calibrations 
between SDSS DR5 and DR7 (R. Tojeiro, private commu- 
nication, see M12 for a detailed discussion). Of the 776,447 
galaxies in our sample, 52,016 galaxies, constituting 6.7 per 
cent of the sample, have VESPA-derived masses only in a 
coarse bin in the time range > 0.42 Gyr (i.e., in bin 29). 
This precludes using these galaxies, as we cannot know how 
to distribute these masses between the 0.42-2.4 Gyr and 
> 2.4 Gyr bins. Furthermore, we exclude those galaxies for 
which the vespa SFH has a reduced \ 2 value in the range 
X 2 > 10, are at redshift z > 1, or that have less than 500 
good pixels in their spectra, according to the MASK array 
included by SDSS for each spectrum [we use only those pix- 
els that have a '0' (good) or '40000000' (emission line) bit 
flag]. These cuts reduce our sample to 707,792 galaxies. 

Of the galaxies in the sample, 29,066 galaxies have more 
than one spectrum, observed on different dates. Of these 
galaxies, after applying the above cuts, we consider in our 
rates analysis only the spectra that are spaced at least 60 
days apart, in order to assure independent search epochs. 
Five of the host galaxies of the SNe la that we present 
below have multiple observations (namely, 0394-51812-554, 
0418-51884-144, 1298-52964-304, 1782-53383-517, and 2019- 
53430-010). Throughout this work we refer to SNe discov- 
ered in SDSS spectra according to the plate, modified Julian 
date (MJD), and fiber in which they are discovered (e.g., 
0814-52443-249 was observed on MJD 52443, or June 18 
2002, in the galaxy targeted by the 249th fiber of plate 814, 
where each SDSS plate contains 640 fibers). The above cut 
leaves 26,592 multiply-observed galaxies: 12,991 with two 
observations; 197 with three epochs; 11 with four; three with 
five; and two galaxies with six epochs. 

We classify the galaxies in our sample according to 
their specific star- formation rate (sSFR), which is com- 
puted as the vespa mass formed in bins 24+25, divided 
by 420 Myr, and divided by the total formed mass. Fol- 
| |2006b[ ), we separate the galaxies 
^ —9.5], star- 
logtsSFR/yr" 1 ) < -9.5], and passive 
< —12] galaxies. As we show in the upper 



lowing Sullivan et al 



into highly star-forming [log(sSFR/yr 
forming [—12 ^ i ... ; ^ i i> , ,. 
[logtsSFR/yr" 1 
panel of Fig. [I] our galaxy sample is dominated by mas- 
sive, old galaxies, and consists of 54.5, 44.7, and 0.8 per 
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Figure 1. Total stellar mass (upper panel) and redshift (lower 
panel) distributions for all galaxies in our sample (solid curve) 
and SN la host galaxies (dashed curve). 



cent passive, star-forming, and highly star-forming galaxies, 
respectively. The SN la host galaxy sample is similarly di- 
vided into 47.8, 52.2, and per cent passive, star-forming, 
and highly star-forming galaxies, respectively. Since only 0.8 
per cent of our galaxies are highly star-forming, we combine 
the star-forming and highly star-forming galaxies into one 
category, which we label 'star-forming galaxies' (i.e., with 
log(sSFR/yr _1 ) > -12). The lower panel of Fig. [I] shows 
the redshift distribution of our galaxy sample and SN la 
host-galaxy sample. The median redshifts of the galaxy and 
SN la host galaxy samples are z = 0.11 and z = 0.09, respec- 
tively. Due to the limits of our detection and classification 
process, we do not find SNe la above z ~ 0.2. 



3 SUPERNOVA DISCOVERY AND 
CLASSIFICATION 

3.1 Method 



1 http:/ /www-wfau. roe. ac.uk/vespa/ 



The galaxy spectra are first de-reddened according to 
Cardelli, Clayton, & Mathisl JT9891 and the E(B - V) val- 
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Table 1. Transient classification criteria 



Transient 


SVD 


SNID 


SNID% 


la 


la 


la 


Ia^ 0.5 


Ia/Ic 


la Ib/c 


la or Ib/c 


Ia> 0.5 


Ib/c 


Ib/c 


Ib/c 


Ib/c> 0.5 


Ic/Ia 


la or Ib/c 


la or Ib/c 


Ia< 0.5 


lib 


Ib/c or lib 


Ib/c or lib 


Ib/cjs 0.5 


II 


IIP, IIL, or Iln 


IIP, IIL, or Iln 


IL> 0.5 


LBV 


LBV 


LBV 




AGN a 


AGN or other 


AGN or other 





Note - In order to be classified into one of the SN subtypes 
in column 1, a candidate must meet both the SVD and 
SNID criteria in columns 2 and 3, and the fraction of best- 
fitting SNID templates in column 4. For example, if a 
candidate is fit as a SN IIP in the SVD stage, and as a SN 
IIL by SNID, and the fraction of best-fitting SNID templates 
is ^ 0.5, that candidate would be classified as a SN II. 
Otherwise, it would be discarded. 

a AGNs were not required to be classified as such by both 
the SVD and SNID stages. 

ues listed in the SDSS spectrum headers. The spectra are 
then masked according to the MASK array included for each 
spectrum. 

In order to identify SNe in the SDSS DR7 galaxy spec- 
tra, we first subtract a model of the galaxy constructed from 
ten SDSS DR2 galaxy eigenspectra ( Yip et al.|2004a l, using 
singular value decomposition (SVD). At this stage, we dis- 
card spectra with best-fitting galaxy models with reduced 
Xr ^ 1, as these spectra either do not host a transient, or 
have a too low signal-to- noise ratio (S/N). The residual spec- 
trum obtained by subtracting the galaxy model is fit with a 
straight line in fx vs. wavelength, which is then subtracted. 
This 'flattened' residual is searched for 'features' by count- 
ing the number of continuous pixels above (or below) zero. 
The residual spectrum must have at least ten features, each 
^ 30 pixels wide. These values were found to be optimal for 
discovering faint SNe, while reducing the number of false 
positives. 

Every galaxy spectrum that satisfies the above criteria 
is then re-analysed by adding SN templates to the SVD- 
constructed galaxy model. We use the sp ectral library from 
the Supernova Identification code (sNirj^J Blondin & Tonry 
|2007[ ), which contains SN (la and core collapse), active galac- 
tic nucleus (AGN), luminous blue variable (LBV), flaring M 
dwarf, and galaxy spectra. For each SVD fit we compute a 
figure of merit, xL composed of the reduced \ 2 value, di- 
vided by the wavelength range covered by the specific tran- 
sient template. This is necessary, as some transient tem- 
plates cover small wavelength ranges due to the instrumental 
setups with which they were taken. Such templates would re- 
ceive a low x 2 value if they fit specific features in the residual 
spectrum, even though a different template from the same 
transient class that covers a larger wavelength range may not 
fit the rest of the residual. This can lead to the misclassifi- 
cation of transients (e.g., SNe la as core-collapse SNe). To 
overcome this problem, we keep only those transient tem- 
plates in the snid library that cover a wavelength range of 
at least 2,900 A, reducing the snid spectral library from 349 
objects with 3,056 spectra to 304 objects with 2,952 spectra. 
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Figure 2. A SN la discovered in the SDSS galaxy spectrum 
0424-51893-355 (plate-MJD-fiber) at z = 0.054. Top: the origi- 
nal galaxy spectrum (black) is fit with a galaxy model (green) 
composed of ten eigenspectra, resulting in a reduced \r = 7.6. 
When fit with both eigenspectra and transient templates (blue), 
a SN la template produces a reduced \r = 1-1- Bottom: residual 
spectrum (black) produced by subtracting the galaxy model from 
the original data. The best-fitting SN la template is shown in red. 



Only transient models with a reduced Xr value smaller than 
the original-galaxy-fit Xr are considered for further analysis. 

The residual spectrum with the best minimal x\ i s m ~ 
put to snid for a second, independent classification, where 
the following criteria must be met: 

(i) The residual must fit more than one template. This elim- 
inates a large number of low-S/N false positives. 

(ii) The best-fitting template must have rlap^ 5, and lapjs 
0.4. The 'lap' parameter quantifies the amount of overlap be- 
tween the residual spectrum and the transient templates at 
the correlation redshift, while 'rlap' is a quality parameter, 



http: / / marwww.in2p3.fr/~blondin /softwarc/snid / index.html 



where higher values denote better correlations (see Blondin 
fc Tonry|2007| for a full description) . 

(iii) The best-fitting templates must be at a redshift within 
Az = 0.015 of the SDSS-measured redshift. 

A large number of false positives are culled at this stage, as 
some residual spectra with low S/N ratios might be fit with 
transient templates using SVD, but will be disqualified by 
snid. 

Finally, the candidate is classified into eight transient 
categories, according to the SVD and snid best-fitting tem- 
plates, and according to the percentage of snid templates 
belonging to a specific category, which we call snid%. The 
transient candidates are classified as la, Ib/c, lib, II, LBV, 
or AGN. We also include two additional categories: Ia/Ic 
(Ic/Ia), which include SNe la (Ic) that were misclassified 
as SNe Ic (la). Were Ia/Ic SNe to be added to the SN la 
sample, they would carry a systematic uncertainty (e.g., if 
two SNe Ia/Ic are discovered, they would add a —2 system- 
atic error to the number of SNe la). However, as detailed in 
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Table 2. SN la population fractions and luminosity functions 



Table 3. Classification purity 



SN la subtype 


fraction 


adopted LF 






[Vega mag] 


normal la 


0.761 


-19.42 ±0.51 


Ia-91bg 


0.024 


-18.30 ± 0.52 


Ia-91T a 


0.197 


-19.90 ± 0.53 


Ia-02cx 6 


0.018 




IIP 


0.394 


-16.41 ±1.23 


IIL 


0.275 


-18.19 ±0.64 


Iln 


0.230 


-17.61 ± 1.61 


lib 


0.101 


-17.40 ± 1.30 


lb 


0.324 


-17.76 ±0.41 


Ic 


0.528 


-16.79 ± 1.28 


peculiar Ib/c 


0.148 


-16.25 ± 1.21 



"This category includes SN2000cx and SN2006gz. 
6 As SN2002cx and SN2005hk were both as faint 
as SN-91bg-like SNe la, they were given its LF. 



Section |4j we find no such SNe. The classification rules are 
summarized in Table [T] 

Along with the SNe in our sample, we also find ~ 980 
AGNs. In order to check whether any of these AGN-hosting 
galaxies also host other types of transients, we re-analyse 
them one more time, using a combination of the galaxy 



eigenspectra, the Vanden Berk et al. (2001 1 composite AGN 



template (which is also included in the snid spectral library), 
and all other transients in the snid spectral library. In other 
words, whereas in the previous stage the data was fit with 
either an AGN or a different transient template, in this stage 
it is fit with both simultaneously. We find no other transients 
in these galaxies. 

We do not mask galaxy emission lines in the first stage 
of analysis, as we have found that the emission lines have 
a negligible effect on our discovery success rate. However, 
since the strength of individual emission lines changes from 
galaxy to galaxy, they are not perfectly fit by the eigen- 
spectra, contaminating the residual spectrum and causing 
snid to misclassify SN candidates. In the re-analysis stage, 
we therefore mask the emission lines of Ha, H/3, [OI] A6300, 
[Oil] A3727, [OIII] A4363, 4959, 5007, Nil A6548, 6583, and 
the Na D A5890, 5896 absorption feature. In Fig. [2] we show 
an example of a SN la discovered by our methodology. 



3.2 Comparison with previous spectral supernova 
searches 

As in this work, M03 and Kll used galaxy eigenspectra, ob- 
tained through principal component analysis, to subtract a 
galaxy model from the data. The residual spectrum, how- 
ever, was analysed differently. M03 performed a wavelet 
transform on the residual to separate the SN signal from 
the noise. This resulted in a 'noiseless' signal which they 
compared, through cross-correlation, to a small set of SN 
la spectral templates. Such manipulation of the data car- 
ries the risk of introducing artificial features, and indeed we 
believe that some of the SNe reported by M03 may have 



been false detections (see below). Kll used the same Yip 
|et al.| ( |2004a| ) eigenspectra as in this work, together with an 
AGN template composed of the first two QSO eigenspectra 



from Yip et al. ( 2004b I, and the Nugent, Kim, & Perlmutter 
( |2002[ > SN la template. Their comparison to SN templates 
was done solely through SVD. 



- 

J: 

5 



la 



Planted as: 
Ib/c 



II 



la 

Ia/Ic 
Ib/c 
Ic/Ia 
II 



99.13 
0.68 
0.09 
0.08 
0.02 



1.4 
5.6 
91.7 
1.2 
0.1 



0.10 
0.24 
0.51 
0.11 
99.04 



Note — All values are percentages. 

Our method differs in four main ways: 

(i) We use the snid as well as the SVD-derived \ 2 criterion. 
This allows us to cull false SN la candidates (as detailed in 
Section [3. 1[ ) and other false positives. 

(ii) We have found that using an AGN template together 
with the galaxy eigenspectra 'robs' some of the continuum 
from a transient in the galaxy spectrum, making it harder 
to discover fainter SNe. Instead, as detailed in Section [3.1[ 
we first find ~ 980 AGN-hosting galaxies, and search them 
for additional transients. By limiting the combination of the 
AGN template and the galaxy eigenspectra to these AGN- 
hosting galaxies alone, we assure that the AGN template 
does not affect the shape of the residual spectra in non- 
AGN-hosting galaxies. 

(iii) By using snid, and the snid templates in the SVD 
phase, we can discover other types of transients, from core- 
collapse SNe to LBVs. Furthermore, snid and the SVD 
phase are modular: new spectra can be added as they be- 
come available (e.g., for new types of transients, such as the 
super-luminous SNe). 

(iv) The data are not manipulated in any way (e.g., there 
is no binning). 

We recover all of the eight examples reported by Kll 
out of the 52 SNe la they discovered, except 0472-51955- 
247 (see their figure 8) which is hosted by a galaxy that 
is not included in our galaxy sample. Of the 19 SNe la re- 
ported in M03, we recover ten. Two of the 19 SNe are hosted 
by galaxies that are not included in our galaxy sample, as 
they do not have associated vespa SFHs. Of the remain- 
ing seven SNe la, the residual spectra of the four galaxies 
hosting SN2000ga, SN2001ki, SN2001kl, and SN2001kr have 
no discernible SN la features, and we deem them false de- 
tections. SN2001kn and SN2001ks may be real SNe la, but 
were not detected by our methodology due to the low S/N 
of their residual spectra. Finally, SN2001kq is detected with 
a redshift that is higher by Az > 0.015 than that of its host 
galaxy, and so is rejected by our detection and classification 
process. The redshift difference appears to be due to a rather 
poor best fit to a noisy residual spectrum, indicating a pe- 
culiar SN la. Differences in the SN la rate analysis between 
this work and M03 and Kll are discussed in Section[5]3] 



3.3 Detection and classification efficiency 

SNe may be missed by our discovery and classification pro- 
cess for numerous reasons, including poor data (such as spec- 
tra with large masked areas), low S/N, low contrast relative 
to the galaxy light, or failure to meet our classification crite- 
ria. In order to quantify these and other systematic effects, 
we simulate the discovery and classification efficiency of our 
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Figure 3. SN la detection efficiency as a function of SN la r J - 
band magnitude. Filled circles denote the fraction of fake SNe la 
discovered in all galaxies in 0.5-mag-wide bins. Error bars indi- 
cate ltr binomial uncertainties. Curves are cubic spline fits to the 
detection-efficiency measurements of SNe la in host galaxies with 
r'-band magnitudes in different ranges, as marked. 
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Figure 4. Distributions of differences between the planted and 
measured </-, r'-, and i'-band magnitudes of the fake SNe la (filled 
squares) and their host galaxies (filled circles), respectively. The 
markers denote the medians of the distributions, and the vertical 
error bars denote the 16th and 84th percentiles. 



method by planting fake SN spectra in a random selection 
of real galaxy spectra. 

The fake SN spectra are drawn from the the subsample 
of snid spectra spanning 2,900 A. In order to simulate the 
real distribution of SNe la in our sample, the host galaxies 
are chosen according to their sSFR. As the rate of SNe la is 
higher in more luminous galaxies, we favor in our simulation 
the selection of host galaxies based on luminosity, according 
to the dependence of the rate of SNe la on host-galaxy r'- 
band luminosity in figure 12 of |Yasuda fc Fukugita| ( |2010[ ). 
Furthermore, as the SN la rate in passive galaxies is ~ 1/3 



the rate in star-forming galaxies (see figure 6 in Sullivan 
|et al.|[2~006b[ ), 25 (75) per cent of the host galaxies in the 
fake sample arc chosen to be passive (star-forming) galax- 
ies. Next, the SN la spectra are planted according to the 
luminosity functions (LFs) and subtype populations from 



Li et al. (2011a I. The Li et al. (2011a I normal-la LF is cen- 



tred on —18.67 mag (Vega), with a standard deviation of 
0.51 mag. This LF is fainter than most measured SN la 



LFs, so we use instead the Yasuda & Fukugita (20101 value 
of —19.42 mag (Vega). We use the difference between these 
values (0.75 mag) to similarly brighten the mean values of 
the |Li et al.] ( |2011a[ ) SN-subtype LFs. Table[2]lists the popu- 
lation fractions and adopted LFs we use for each SN subtype. 
Varying the mean brightness of the LFs (e.g., using the orig- 
inal Li et al.|2011a values) does not affect the overall shape 
of the resulting detection efficiency function. 

For each fake spectrum we choose a random peak B- 
band magnitude from the appropriate LF, and translate it 
into a stretch value, s, according to the |Perlmutter et aL] 



(1999) stretch-luminosity relation 



nis.s = rriB,\ + a(s — 1), 



(1) 



where m_g, s is the peak B-band magnitude of the stretched 
light curve, ms,i is the peak _B-band magnitude of a normal 
s = 1 SN la, and a = 1.52 dAstier et al.||2006h. From this 



stretched light curve we choose a random epoch, and assign 
the fake spectrum the r'-band magnitude on that epoch. 



Following Sullivan et al. (2006a I, we limit the stretch range 



to 0.6 < s < 1.4. 

We conduct three separate simulations, in each of which 
~ 10,000 fake spectra of a specific SN subtype are planted 
in ~ 75,000 random galaxy spectra. The SN la simu- 
lation includes all normal SNe l a, SN-91T-like S Ne [in- 
cluding SN2 000CX ( |Li et al.||200"Tj ) and SN2006gz ( |Hicken| 



et al. 20071], SN-91bg-like SNe, and SN2002cx-like SNe 



SN20 02CX flLi et al.||2003| ) and SN2005hk ( |Phillips 
et al.||2007|>]. We do not include t he peculiar SNe SN2 002ic 
( [Haniuy et al"1|2003[), SN2003 fg dHowell et al]|2006 |), and 



SN2005gj ( |Prieto et aLp 007). The SN Ib/c simulation in- 
cludes normal SNe lb, SNe Ic, and SNe lib, along with pe- 
culiar SNe lb and broad-lined SNe Ic. The SN II simulation 
includes all SN IIP, SN IIL, and SN Iln spectra included in 

SNID. 

Similarly to the SNe la, the SN Ib/c and SN II fakes 
are planted according to the LFs and population fractions 



presented in tables 6 and 7 of Li et al. ( 201 la I. We apply the 



same magnitude correction of 0.75 mag to these LFs, which 
are presented in Table [2] Contrary to the fake SNe la, the 
fake SNe Ib/c and SNe II are only planted in star-forming 
galaxies, i.e., those galaxies with log^SFR/yr" 1 ) ^ —12. 

Table[3]lists the purity of our classification method, i.e., 
what fraction of fake SNe were classified into the different 
subtypes listed in Table [T] The table shows that our code is 
efficient at classifying SNe la and SNe II. SNe Ib/c, however, 
are harder to classify, as they may appear similar to SNe la 
(especially SNe Ic that exhibit no hydrogen or helium in 
their spectra). There is a 1.4-per-cent chance that one of 
our SNe la is in fact a misclassified SN Ib/c. We take this 
into account by adding a systematic uncertainty of — 1 to our 
sample of 90 SNe la, described in Section [4] Furthermore, 
there is a 5.6-per-cent probability that a SN Ib/c will be 
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Figure 5. Left and centre: distributions of differences between the planted and recovered SVD- and SNlD-derived ages of the fake SNe la, 
respectively, in 1-day-sizcd bins. Right: comparison of the SVD- and SNlD-derived ages (crosses). The solid curve represents the 1:1 line. 



misclassified as a Ia/Ic. However, we have no such SNe in 
our SN la sample. 

We present the results of our SN la simulation in Fig- 
ures [3Jj5] Fig. [3] shows our SN la detection efficiency as a 
function of the r'-band apparent magnitude of the planted 
SNe. We reach 50 per cent efficiency at r' — 20.7 mag. Since 
the brightness of the host galaxy affects the probability of 
detecting the SN (i.e., it is easier to discover SNe in fainter 
galaxies, and vice versa), we divide the SN la detection effi- 
ciency measurements shown in Fig. [3] into four subsets, ac- 
cording to the host-galaxy r'-band magnitudes. These sub- 
sets are: r> h > 19, 18 < r' h sC 19, 17 < r\ < 18, and r' h sC 17 
mag. To the efficiency curves of each of these subsets we fit 
a cubic spline. 

In Fig. [4] we show the distributions of the differences be- 
tween the planted and recovered apparent magnitudes, Am, 
of the fake SNe la, and their host galaxies, in the g' , r 1 , and 
i bands. Because the transient spectra are not orthogonal 
with the set of galaxy eigenspectra, there is a flow of flux 
between them, which causes the SNe to appear brighter, on 
average, in the g' band, and fainter in the r' and i' bands. 
The host galaxies, conversely, appear fainter and brighter, 
respectively, in these bands, as shown. Throughout this pa- 
per we report SN and host-galaxy apparent magnitudes in 
the r' band, where the dispersion in Am is smallest. More- 
over, as the average offset in the measured r'-band magni- 
tude of the SNe la caused by the flux flow is negligible, and as 
the uncertainties in brightness are smaller than the 0.5- mag- 
wide bins we use for the detection efficiency measurements 
of Fig. |3j we do not consider this offset any further. 

The left and centre panels of Fig. [5] show the distri- 
butions of differences between the planted and SVD- and 
SNiD-derived ages, respectively, in 1-day-sized bins. The right 
panel of Fig. [5] shows a comparison of the SVD- and snid- 
derived ages. While the derived ages are almost always 
within ±1 day of the planted values, the comparison be- 
tween the SVD- and SNiD-derived ages in the right panel 
of Fig. [5] shows that there is a larger source of systematic 
uncertainty contributing to the discrepancy between them. 
If we assume that the overall uncertainty is identical for all 
the SNe la, then we can estimate it by fitting the SVD- 



to the SNiD-derived ages and setting x = 1- After limiting 
the ages of the fake SNe la to the range of the real SNe la 
in our sample, —10-55 days relative to _B-band maximum 
light, we estimate the overall uncertainty in the age of the 
SNe la to be ±6 days. In a similar manner, we estimate the 
uncertainty in the ages of the SNe IIP to be ±33 days. 



3.4 Supernova age and stretch recovery efficiency 

In order to measure the stretches of the SNe la in our sam- 
ple, we substit ute the snip tem plates used in the discovery 
phase with the Guy et al. (20071 salt^] templates. Similarly 



to the Yip et al. (2004a I eigenspectra, the salt2 templates 
are composed of two base spectral sequences: Mo(p, A) and 
Mi (p, A) , where p is time since maximum light in the rest- 
frame B-band, A is the rest-frame wavelength, Mo(p, A) is 
the average time sequence of eigenspectra, and Mi (p, A) is 
the spectral sequence of variations induced by the stretch of 
the SNe. Following equation 1 in 
the SALT2 eigenspectra linearly: 



Guy et al. (20071, we add 



F(p, A) = x x [M (p, A) + xiMi(p, A)], 



(2) 



but without the exponential colour correction. The variable 
xi can be converted into the stretch parameter s by using 



the Guy et al. (20071 transformation 



s = 0.98 + 0.091a;! + 0.003x1 - 0.00075x?. (3) 

We first gauge the accuracy of our stretch measurements 
by planting fake SNe la made up of the salt2 templates with 
predetermined stretches. Fig.[6]shows the distribution of dif- 
ferences between the planted and recovered stretch values for 
a sample of 7,200 fake SNe la planted using the salt2 tem- 
plates, in bins of As = 0.1. The distribution is more sharply 
peaked than a Gaussian, and the range of stretches enclos- 
ing 68 per cent of the area is ±0.03. We take this to be the 
uncertainty of our stretch measurements. As in Section [3. 3[ 
fitting the SALT2-derived ages to those derived by SVD re- 
sults in the same uncertainty of ±6 days. 



http: / / supernovae.in2p3.fr / ~guy / salt / 
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Figure 6. Results of using our discovery code with the 
SALT2 templates: the distribution of differences between the 
planted and recovered stretch values for a sample of 7,200 fake 
SNe la planted using the SALT2 templates. 
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Figure 7. Stretch values for 90 SNe la, as measured using the 
SALT2 templates. The stretch values in the upper (lower) panel 
are for the 47 (43) SNe la in star-forming (passive) galaxies, where 
star-forming (passive) galaxies are chosen to have log(sSFR/yr _1 ) 
above (below) — 12. 



4 SUPERNOVA SAMPLE 

We find a total of 100 SNe, of which 90 are SNe la and 10 
are SNe IIP. Table [| lists the SNe and their properties. The 
residual spectra and best-fitting SN templates are shown in 
Figures [8] and [9] Each SN is listed according to the plate, 
MJD, and fiber in which it was discovered. We measure the 
r'-band magnitudes of the SNe by performing synthetic pho- 
tometry on the SN residual spectrum, as obtained by sub- 
tracting the galaxy model from the original data. Similarly, 
the r'-band magnitudes of the SN host galaxies are mea- 
sured using synthetic photometry on the original data, after 
subtracting the SN residual. Finally, we measure the stretch 
of each SN using the salt2 templates, as detailed in Sec- 
tion E! 

Apart from these SNe, we detect ~ 980 AGNs and 
several dozen transients which might be either SNe Iln or 
AGNs. Some SN Iln spectra, which exhibit broad hydrogen 
lines, are indistinguishable, without further spectroscopic or 
photometric followup, from AGNs, which have similar emis- 
sion lines. 

Some of the SNe la in our sample are sufficiently bright 
so as to significantly distort the shape of the combined 
galaxy+SN spectrum. This results in a systematically bi- 
ased vbspa fit and an incorrect SFH for the host galaxy. 
In order to correct this, we replace the SFHs of such host 
galaxies with those of 'surrogate' galaxies that did not host 
SNe, by fitting the continuum of the SN host galaxy with 
those of all galaxies within our galaxy sample that have the 
same redshift, within Az — 5 x 10~ 4 , and r'-band magni- 
tude within Am = 0.05 mag. We inspect by eye the resulting 
five best-fitting galaxies. In most cases, we select the galaxy 
with the minimal reduced Xr value as the surrogate galaxy. 
In some cases, a better fit, with a Xr value that is slightly 
larger (up to two tenths) of the minimal Xr> ls discernible by 
eye, and we choose it instead. We find that no surrogates are 
necessary for host galaxies where the contrast between the 
host and the SN la is ^ 0.1. For SN-hosting galaxies with 
multiple, SN-free observations (i.e., 0394-51812-554, 0418- 



51884-144, 1298-52964-304, and 2019-53430-010), we use the 
formed masses derived from the vbspa fits to the SN-free 
galaxy spectra. 

Figure[7]shows the stretch distributions of the 90 SNe la 
in our sample. Many previous studies have found that star- 
forming galaxies tend to host higher-stretch SNe la than 
passive galaxies (e .g., |Hamuy et al.||2000| |Neill et al.||2009| 



Hicken et al.|2009[ ). For example, |Sullivan et al.| ( |2006b[ ) find 
a formal difference between the stretches of SNe la hosted 
by passive and star-forming galaxies at z > 0.1, with stretch 
medians of ~ 0.93 and ~ 0.98, respectively. Our sample has 
stretch medians of 0.95 (1.0) for the 43 (47) SNe la hosted by 
passive (star-forming) galaxies, which is a stretch difference 
of the same sense and similar magnitude to those seen by 
previous studies. However, a number of statistical tests (x 2 , 
Kolmogorov-Smirnov, Student) all indicate that, in the case 
of our sample, these differences are not significant. Whereas 
the apparent statistical accuracy of our stretch measure- 
ment, ±0.03, is similar to that of Sullivan et al. (2006b I, 



±0.01-0.02, our measurements might suffer from systematic 
or larger statistical uncertainties currently unaccounted for. 

When using the figure of merit x\ introduced in Sec- 
tion |3T] we optimize our detection and classification method 
towards the discovery of SNe la, and recover only two of the 
ten SNe IIP in Table [1J due to the relatively smaller wave- 
length ranges covered by the snid SN IIP templates. This, 
together with the overall small size of the SN IIP sample, 
precludes deriving a SN IIP rate in the context of this work. 



5 SUPERNOVA RATE ANALYSIS 
5.1 Visibility time 

In order to measure SN la rates, we calculate a visibility time 
for each galaxy in our sample. This is the period of time dur- 
ing which we could have detected a SN la at that galaxy's 



Supernovae in SDSS DR7 Spectra 



Table 4. SNe discovered in SDSS DR7 



SDSS ID 


Date 


a 


S 


z 


Agci 


Age 2 


s 


r'sN 


r'n 


Xgal 


*SN 


Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


0271-51883-171 


05/12/00 


154.5020 


-0.0328 


0.065 


27 


42 


0.91 


19.54 


18.11 


2.4 


1.2 


la 


0291-51928-076 


19/01/01 


191.8892 


0.0992 


0.086 


28 


21 


1.06 


20.03 


17.87 


1.7 


1.0 


la 


0313-51673-154 


09/05/00 


231.8958 


-0.0595 


0.044 


33 


38 


1.11 


19.85 


17.43 


1.6 


1.1 


la 


0328-52282-570 


08/01/02 


175.9218 


-1.4770 


0.125 


27 


30 


0.84 


19.94 


17.99 


1.5 


1.1 


la 


0358-51818-181 


01/10/00 


263.1189 


56.0737 


0.123 


21 


30 


1.05 


20.19 


18.21 


1.7 


1.1 


la 


0394-51812-554 


25/09/00 


14.0069 


0.4353 


0.146 


-5 


-5 


1.12 


20.03 


18.08 


1.3 


1.0 


la 


0418-51884-144 


06/12/00 


9.8403 


14.4699 


0.017 


12 


16 


1.04 


19.46 


17.36 


2.1 


1.1 


la 


0424-51893-355 


15/12/00 


19.6493 


14.6835 


0.054 


7 


7 


0.94 


17.97 


16.98 


7.6 


1.1 


la 


0438-51884-166 


06/12/00 


121.0888 


46.7870 


0.187 


16 


12 


0.88 


21.30 


18.16 


1.1 


1.0 


la 


0438-51884-462 


06/12/00 


120.8026 


47.6138 


0.117 


2 


-5 


0.91 


19.54 


18.06 


2.2 


1.0 


la 


0452-51911-319 


02/01/01 


140.6214 


57.9081 


0.063 


9 


10 


0.83 


19.43 


17.93 


1.8 


1.1 


la 


0480-51989-024 


21/03/01 


147.9711 


1.1016 


0.063 


2 


-8 


1.00 


18.77 


17.77 


4.4 


1.5 


la 


0498-51984-102 


16/03/01 


212.7430 


64.8475 


0.141 


20 


18 


0.99 


20.43 


18.09 


1.2 


0.9 


la 


0500-51994-100 


26/03/01 


149.8156 


0.9673 


0.088 


11 


9 


0.98 


19.64 


18.43 


1.8 


1.0 


la 


0578-52339-314 


06/03/02 


159.7062 


4.0156 


0.129 


6 


3 


1.04 


19.54 


17.87 


1.8 


1.2 


la 


0604-52079-209 


19/06/01 


204.9366 


62.3936 


0.136 


2 


-5 


0.95 


20.11 


18.38 


1.2 


1.0 


la 


0606-52365-412 


01/04/02 


213.3167 


62.0323 


0.142 


6 


-2 


0.88 


20.15 


17.45 


1.4 


1.2 


la 


0622-52054-011 


25/05/01 


244.3058 


48.4744 


0.104 


16 


17 


1.07 


19.88 


18.21 


1.7 


1.3 


la 


0738-52521-360 


04/09/02 


336.7925 


13.6661 


0.150 


27 


22 


1.04 


20.19 


18.14 


1.1 


1.0 


la 


0745-52258-092 


15/12/01 


350.7757 


13.6668 


0.041 


54 


68 


0.89 


20.16 


17.40 


1.1 


0.9 


la 


0762-52232-067 


19/11/01 


129.7765 


43.7123 


0.125 


31 


34 


1.06 


20.39 


18.03 


1.0 


0.9 


la 


0814-52443-249 


18/06/02 


243.4005 


43.7256 


0.112 


8 


7 


0.96 


20.01 


18.18 


1.5 


0.9 


la 


0844-52378-462 


14/04/02 


184.4150 


5.3234 


0.104 


31 


36 


1.05 


20.32 


17.85 


1.2 


1.0 


la 


0905-52643-213 


04/01/03 


157.4940 


53.5018 


0.137 


42 


38 


1.02 


20.60 


18.25 


1.0 


0.9 


la 


0915-52443-543 


18/06/02 


211.3174 


-1.7113 


0.054 


29 


24 


1.03 


19.33 


17.44 


1.5 


0.9 


la 


0966-52642-221 


03/01/03 


172.2523 


48.7331 


0.074 


3 


2 


0.89 


19.46 


18.26 


1.8 


1.1 


la 


1038-52673-135 


03/02/03 


191.8526 


53.7308 


0.153 


12 


11 


0.88 


20.59 


19.08 


1.2 


1.0 


la 


1059-52618-144 


10/12/02 


116.8937 


27.4465 


0.063 


33 


32 


1.04 


19.43 


17.51 


1.9 


0.9 


la 


1171-52753-185 


24/04/03 


244.8403 


41.0899 


0.038 


10 


8 


0.89 


18.00 


16.16 


3.9 


0.9 


la 


1189-52668-239 


29/01/03 


132.2645 


5.8377 


0.126 


3 


4 


0.78 


20.01 


18.07 


1.4 


1.0 


la 


1205-52670-632 


31/01/03 


122.8267 


26.1661 


0.144 


12 


13 


0.97 


20.25 


17.58 


1.3 


1.1 


la 


1278-52735-425 


06/04/03 


189.1061 


50.6116 


0.106 


6 


6 


0.95 


19.51 


17.64 


1.5 


1.0 


la 


1289-52734-413 


05/04/03 


220.2447 


45.1308 


0.074 


8 


8 


0.93 


19.53 


16.57 


1.0 


0.8 


la 


1298-52964-304 


21/11/03 


129.7903 


7.4088 


0.046 


32 


48 


0.96 


18.92 


17.47 


2.0 


1.0 


la 


1310-53033-459 


29/01/04 


173.5529 


58.2622 


0.122 


16 


14 


1.06 


19.56 


17.67 


1.9 


0.9 


la 


1324-53088-169 


24/03/04 


211.9955 


54.3631 


0.067 


52 


52 


1.06 


19.87 


17.60 


1.5 


1.0 


la 


1337-52767-480 


08/05/03 


245.0462 


38.1115 


0.130 


2 


4 


0.84 


19.80 


18.12 


1.5 


1.0 


la 


1392-52822-147 


02/07/03 


240.3856 


26.9875 


0.050 


14 


14 


0.93 


17.84 


16.90 


6.4 


1.1 


la 


1400-53470-234 


10/04/05 


230.1879 


36.8118 


0.103 


6 


3 


1.03 


18.89 


16.58 


1.8 


0.9 


la 


1400-53470-351 


10/04/05 


229.2268 


37.1240 


0.116 


-9 


-14 


1.16 


19.66 


17.21 


1.1 


0.9 


la 


1403-53227-456 


10/08/04 


237.2391 


33.9571 


0.128 


27 


33 


0.84 


20.26 


18.39 


1.2 


1.0 


la 


1445-53062-067 


27/02/04 


177.1144 


42.1320 


0.086 


3 


-3 


0.89 


19.65 


17.41 


1.3 


0.9 


la 


1462-53112-638 


17/04/04 


202.1417 


41.8523 


0.028 


8 


7 


0.99 


16.52 


16.53 


13.3 


1.1 


la 


1581-53149-470 


24/05/04 


235.1032 


32.8659 


0.054 


30 


40 


0.90 


18.97 


17.69 


3.3 


1.0 


la 


1598-53033-380 


29/01/04 


155.7088 


11.7030 


0.102 


2 


-1 


0.92 


19.67 


17.99 


1.6 


1.0 


la 


1645-53172-349 


16/06/04 


217.2314 


35.0823 


0.121 


21 


34 


1.14 


19.97 


16.95 


1.3 


1.0 


la 


1697-53142-506 


17/05/04 


199.1256 


12.6769 


0.151 


3 


4 


1.06 


19.92 


17.66 


1.3 


0.9 


la 


1700-53502-302 


12/05/05 


203.6364 


11.1324 


0.095 


3 


1 


1.21 


19.59 


16.96 


1.3 


1.1 


la 


1710-53504-488 


14/05/05 


219.0440 


12.4451 


0.085 


27 


24 


0.90 


19.70 


18.03 


1.2 


1.0 


la 


1744-53055-210 


20/02/04 


149.9506 


11.4737 


0.060 


9 


9 


0.94 


18.37 


19.20 


6.6 


1.1 


la 



(1) - SDSS DR7 plate, MJD, and fiber in which the SN was discovered. 

(2) - Date on which the SN was captured, in dd/mm/yy. 

(3) -(4) - Right ascensions and declinations (J2000). 

(5) - SN host-galaxy redshift. 

(6) -(7) - SN age, in days, as measured by SVD and SNID, respectively. SN la (IIP) ages have an uncertainty of ±6 (33) days. 

(8) - SN stretch, as measured with the SALT2 templates. All stretches have an uncertainty of ±0.03. 

(9) -(10) - SN and host-galaxy r'-band magnitudes. 

(11)-(12) - Reduced x 2 value of galaxy and galaxy+transient fits. 
(13) - SN type. 
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Table 4. SNe discovered in SDSS DR7 - continued 



SDSS ID 


Date 


a 


S 


z 


Agci 


AgC2 




s 






Xgal 

(11) 




Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(12) 


(13) 


1755-53386-309 


16/01/05 


173.1052 


14.6202 


0.082 


21 


21 


1 

± 


1 o 


19.32 


17.13 


1.7 


0.9 


la 


1758-53084-523 


20/03/04 


127.3894 


8.8682 


0.112 


2 


-3 


u 


oO 


19.62 


17.75 


1.6 


1.0 


la 


1782-53383-517 


13/01/05 


126.1666 


54.6689 


0.063 


55 


43 


i 


uz 


20.02 


17.33 


1.5 


1.0 


la 


1788-54468-126 


03/01/08 


145.0818 


63.2817 


0.120 


2 


-1 


i 


nn 
uu 


20.46 


18.02 


1.2 


1.0 


la 


1791-54266-583 


15/06/07 


194.1119 


10.2277 


0.107 


12 


10 


i 


U± 


19.31 


17.65 


2.0 


1.2 


la 


1793-53883-040 


28/05/06 


196.4333 


9.5797 


0.055 


52 


50 


i 


Ul 


19.68 


17.03 


1.4 


0.9 


la 


1801-54156-371 


25/02/07 


202.8033 


7.9573 


0.123 


12 


34 


i 


nn 
uu 


19.81 


17.79 


1.4 


1.1 


la 


1803-54152-260 


21/02/07 


205.4068 


5.8710 


0.060 


33 


40 


n 
u 


OO 


19.80 


17.11 


1.3 


1.0 


la 


1843-53816-491 


22/03/06 


223.2384 


31.0742 


0.094 


55 


50 


i 


Uo 


20.01 


17.73 


1.5 


1.0 


la 


1944-53385-434 


15/01/05 


144.4477 


28.2876 


0.153 


-1 


-5 


n 

u 


OO 


20.10 


18.49 


1.6 


1.0 


la 


1949-53433-080 


04/03/05 


150.8604 


32.1330 


0.166 


3 


-5 


1 


uo 


20.23 


18.08 


1.1 


0.9 


la 


1957-53415-232 


14/02/05 


155.6271 


35.6764 


0.128 


21 


17 


u 


yo 


20.23 


18.00 


1.2 


0.9 


la 


2019-53430-010 


01/03/05 


161.1689 


30.6343 


0.072 


2 


3 


1 

1 


Ul 


18.87 


18.68 


4.8 


1.6 


la 


2159-54328-161 


16/08/07 


230.6900 


19.7057 


0.109 


12 


12 


u 


on 

o\) 


19.70 


17.65 


1.2 


0.9 


la 


2165-53917-406 


01/07/06 


234.5504 


25.0456 


0.067 


27 


35 


1 


07 
U i 


19.22 


17.69 


1.3 


0.9 


la 


2173-53874-154 


19/05/06 


240.4273 


20.5769 


0.123 


2 


-1 


U 


07 
O t 


20.30 


17.80 


1.4 


1.1 


la 


2199-53556-232 


05/07/05 


240.3188 


17.7678 


0.045 


13 


15 


u 


QA 


17.84 


16.65 


10.9 


1.5 


la 


2202-53566-403 


15/07/05 


246.0988 


15.6767 


0.084 


4 


4 


u 


QQ 

yy 


18.54 


18.45 


8.1 


1.2 


la 


2218-53816-295 


22/03/06 


171.5008 


26.0536 


0.158 


6 


9 


1 


1 
Ul 


20.24 


17.96 


1.4 


1.2 


la 


2222-53799-480 


05/03/06 


176.1602 


29.8899 


0.076 


12 


15 


u 


QD 

yu 


19.29 


18.23 


3.3 


1.1 


la 


2376-53770-183 


04/02/06 


158.4883 


20.3405 


0.087 


2 


-3 


1 

1 


09 

uz 


18.88 


16.90 


2.0 


1.0 


la 


2420-54086-142 


17/12/06 


123.7092 


11.7155 


0.088 


-3 


-1 


u 


y4 


20.77 


18.05 


1.5 


1.2 


la 


2429-53799-033 


05/03/06 


132.1170 


14.4232 


0.069 


38 


36 


1 


O 1 } 
Uo 


19.75 


17.71 


1.5 


1.0 


la 


2430-53815-267 


21/03/06 


132.4331 


12.2987 


0.050 


33 


36 


1 


Uo 


19.47 


18.77 


2.4 


1.1 


la 


2499-54176-550 


17/03/07 


172.2386 


17.2243 


0.143 


-3 


-9 


u 


QO 

yu 


20.11 


18.10 


1.3 


0.9 


la 


2594-54177-348 


18/03/07 


158.4080 


16.6030 


0.052 


21 


23 


1 

1 


1 Q 


20.69 


17.48 


1.5 


1.1 


la 


2603-54479-486 


14/01/08 


196.4791 


17.8272 


0.078 


28 


26 


1 


UD 


20.91 


18.56 


1.2 


0.9 


la 


2609-54476-295 


11/01/08 


181.4544 


18.9829 


0.168 


3 


-5 


1 


1 9 
1Z 


19.96 


18.08 


1.6 


1.1 


la 


2614-54481-257 


16/01/08 


190.8844 


18.9626 


0.168 


8 


4 


1 

1 


OO 
UU 


20.42 


18.37 


1.2 


1.0 


la 


2664-54524-468 


28/02/08 


200.7558 


24.5066 


0.073 


6 


6 


n 
u 


yo 


19.11 


17.63 


2.4 


1.2 


la 


2744-54272-561 


21/06/07 


211.5930 


16.4834 


0.014 


40 


45 


ci 
u 


QQ 
yy 


17.40 


15.23 


4.5 


0.8 


la 


2747-54233-613 


13/05/07 


218.8254 


15.5351 


0.107 


12 


10 


u 


ST 

O I 


19.90 


17.88 


1.6 


1.1 


la 


2754-54240-593 


20/05/07 


234.2158 


12.1352 


0.094 


27 


26 


u 


OO 


19.90 


18.09 


1.9 


1.1 


la 


2758-54523-082 


27/02/08 


213.1021 


17.0065 


0.174 


6 


6 


u 


yo 


20.25 


18.28 


1.4 


1.1 


la 


2768-54265-233 


14/06/07 


233.2249 


13.1177 


0.073 


27 


24 


1 


1 
1U 


19.54 


17.29 


1.3 


1.0 


la 


2771-54527-05 


02/03/08 


212.4780 


19.6616 


0.077 


16 


17 


1 

1 


OO 

uy 


20.21 


16.91 


1.2 


1.0 


la 


2792-54556-210 


31/03/08 


226.3821 


17.9847 


0.036 


13 


15 





93 


17.35 


16.58 


15.7 


1.1 


la 


2886-54498-598 


02/02/08 


164.2943 


9.4009 


0.088 


6 


10 





96 


19.86 


18.25 


2.0 


1.3 


la 


2949-54557-440 


01/04/08 


249.2445 


27.1098 


0.140 


21 


22 





89 


20.45 


18.52 


1.2 


1.0 


la 


2954-54561-572 


05/04/08 


233.5161 


2.2105 


0.034 


-1 


5 





83 


18.82 


16.76 


1.3 


0.8 


la 


0437-51869-328 


21/11/00 


119.5556 


44.0190 


0.047 


19 









20.00 


17.48 


1.5 


1.1 


IIP 


0864-52320-082 


15/02/02 


129.9163 


35.6542 


0.160 


4 


5 






20.00 


18.31 


1.2 


1.0 


IIP 


1207-52672-512 


02/02/03 


126.2081 


29.6123 


0.040 


87 


36 






20.23 


18.06 


1.4 


1.1 


IIP 


1231-52725-553 


27/03/03 


186.9265 


9.9579 


0.070 


34 


49 






20.15 


18.41 


1.2 


1.0 


IIP 


1406-52876-528 


25/08/03 


243.2197 


30.8496 


0.048 


34 


49 






20.25 


18.18 


1.3 


1.0 


IIP 


1459-53117-022 


22/04/04 


198.2796 


46.0984 


0.030 


35 


50 






19.88 


19.33 


2.4 


1.7 


IIP 


1684-53239-484 


22/08/04 


245.6866 


32.6592 


0.041 


62 


36 






19.85 


17.01 


1.9 


1.4 


IIP 


1755-53386-516 


16/01/05 


174.8064 


15.0377 


0.014 


63 


82 






19.51 


18.71 


3.2 


1.2 


IIP 


2103-53467-081 


07/04/05 


180.8496 


35.3258 


0.028 


19 


34 






19.22 


17.48 


1.8 


0.9 


IIP 


2593-54175-334 


16/03/07 


157.8944 


19.0687 


0.034 


35 


113 






19.98 


19.68 


2.0 


1.4 


IIP 



(1) - SDSS DR7 plate, MJD, and fiber in which the SN was discovered. 

(2) - Date on which the SN was captured, in dd/mm/yy. 

(3) -(4) - Right ascensions and declinations (J2000). 

(5) - SN host-galaxy rcdshift. 

(6) -(7) - SN age, in days, as measured by SVD and SNID, respectively. SN la (IIP) ages have an uncertainty of ±6 (33) days. 

(8) - SN stretch, as measured with the SALT2 templates. All stretches have an uncertainty of ±0.03. 

(9) -(10) - SN and host-galaxy r'-band magnitudes. 

(11)-(12) - Reduced x 2 value of galaxy and galaxy+transicnt fits. 
(13) - SN type. 
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0271-51883-171: x 2 (Gal)=2.4, % 2 (Gal+SN)=1.2, z=0.065, s=0.91 



0291-51928-076: x 2 (Gal)=1.7, x 2 (Gal+SN)=l, z=0.086, s=l.l 
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0418-51884-144: x 2 (Gal)=2.1, x 2 (Gal+SN)=l.l, z=0.017, s=l 
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0438-51884-166: x 2 (Gal)=l.l, x 2 (Gal+SN)=0.97, z=0.19, s=0.8 
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0328-52282-570: % (Gal)=1.5, jf (Gal+SN)=1.1, z=0.12, s=0,i 
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-0.4 




4000 5000 6000 7000 8000 9000 

Wavelength [A] 

0394-51812-554: x 2 (Gal)=1.3, x 2 (Gal+SN)=l, z=0.15, s=l.l 
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0424-51893-355: x 2 (Gal)=7.6, x 2 (Gal+SN)=l.l, z=0.054, s=0.94 
5 
4 

= 3 
3 2 
1 

4000 4500 5000 5500 6000 6500 7000 7500 
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0438-51884-462: x 2 (Gal)=2.2, x 2 (Gal+SN)=l, z=0.12, 5=0.91 
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Figure 8. SNe la discovered in SDSS DR7 galaxy spectra. The residual spectrum, obtained by first fitting galaxy eigenspectra and 
transient templates to the original spectrum, and then subtracting the resulting galaxy model, is shown in black. The best-fitting SN la 
template is overlaid in red. The flux is in units of 10 -16 erg cm -2 s -1 A" 1 . The title of each panel details the plate, MJD, and fiber 
in which it was discovered; the reduced x 2 value obtained when fitting only galaxy eigenspectra to the spectrum, XrC-»al); the reduced 
X 2 obtained from the best-fitting combination of galaxy eigenspectra and transient templates, xjKGal+SN); the redshift of the SN-host 
galaxy, z; and the stretch parameter, s, obtained by using the SALT2 templates. Similar figures for the rest of our SN la sample are 
available in the electronic version of the paper - see Supporting Information. 
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0437-51869-328: % 2 (Gal)=1.5, x 2 (Gal+SN)=l.l, z=0.047 



0864-52320-082: x 2 (Gal)=1.2, x 2 (Gal+SN)=0.98, z=0.16 
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Figure 9. Same as Fig. [8] for SNe IIP discovered in SDSS DR7 galaxy spectra. 



redshift, with specific stretch and host-galaxy extinction val- 
ues, and the SN la detection efficiency curve we measured 
in Section |3.3| We construct SN la light curves using the 



Hsiao et al. 



1 2007[ ) SN la template spectra. Each spectrum 



is reddened using the Cardelli et al. (1989) extinction law 



and Av values drawn from the Neill et al. ( 2006 1 extinction 



model: a one-sided positive Gaussian centred on Av = 



with a standard deviation of a — 0.62. Next, the spectra 
are redshifted according to the redshift of the galaxy. We 
apply synthetic photometry to the reddened and redshifted 
spectra to construct an r'-band light curve. This light curve 
is then stretched according to Equation [I] Each galaxy in 
our final sample is assigned a random stretch value from a 
Gaussian distribution with a standard deviation of a s = 0.1, 
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and the mean of the Gaussian is selected according to the 
sSFR of the galaxy: fj, s = 1 (fi s = 0.9) for star-forming 
(passive) galaxies ( Sullivan et al.||2006a[). Following Sulli- 
Ivan et al. (2006a I and Meyers et al.l J2012I), the stretch val- 



ues for star-forming (passive) galaxies are restricted to the 
range 0.6 < s < 1.4 (0.6 < s < 1.1). We use the |Yasuda~fc 
Fukugita (20101 LF to assign absolute B-band magnitudes 
at maximum light to each light curve. Each point in the 
light curve is then multiplied by the appropriate detection 
efficiency, according to the r'-band magnitude of its host 
galaxy, as shown in Fig. [3] The visibility time is 



it] 



e{m(t)]dt, 



(4) 



where m(t) is the redshifted and reddened template SN la 
light curve in the r' band, and e(m) is the detection effi- 
ciency, as a function of magnitude. 



5.2 The rate-mass relation 

The mass-normalized SN la rate is the number of SNe la, 
Ni a>i , in galaxies within a specific mass range i, divided by 
the sum of the visibility times, t v ,j, of the n galaxies within 
that mass range, weighted by their monitored stellar mass, 
M*j (i.e., the mass of the galaxy within the spectral aper- 
ture, as observed at its present redshift, after mass loss dur- 
ing stellar evolution): 



Rla 



3=1 



(•») 



For the total stellar mass of the galaxies in our sample (used 
only for the purpose of binning the galaxies into various 
total-mass ranges), we use the vespa values from T09, with 
the 0.55 correction, and with the original T09 aperture cor- 
rection. As in other such studies, we present our rates in 
units of SNuM (10" 12 Mo" 1 yr" 1 ). We measure the mass- 
normalized SN la rates in four total mass bins, where the 
limits of the bins are chosen so that each bin contains ap- 
proximately the same number of SNe la. We measure the 
rates of all galaxies, and specifically of star-forming and pas- 
sive galaxies, distinguished according to their sSFR, using 
the previously defined borders. Over the redshift range of 
our SN la sample, 0.014 < z < 0.19, we do not find a depen- 
dence of the rate on redshift, whether for all galaxies or for 
galaxies separated by sSFR. We therefore analyse together 
all redshifts in our sample. The cited uncertainties of the 
rates are the Poisson uncertainties on the number of SNe la 
in each mass bin. These rates are summarized in Table [5] 

The mass-normalized SN la rate, averaged over 
all masses and redshifts in our sample, is 7?i a ,M = 
0.10 ± 0.01 SNuM. For passive galaxies alone, the mass- 
normalized SN la rate is i?i a , M = 0.082±g:gi| SNuM. 
This last value is consistent with the values ob- 
tained by the two other main ga laxy-targeted surveys , 
0.044±q J&t SNuM for E/S0 galaxies ( jMannucci et al.|2005| , 
and 0.125j:g;^!(0.028) [0.104l£°^(0.024)J HNuM for ellip- 
tical (SO) galaxies (Lll), where the uncertainties are statis- 
tical and systematic, respectively. 

Fig.|10|shows the SN la rate per unit mass, as a function 
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0.2 
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0.05 




All galaxies 
Star-forming galaxie 
Passive galaxies 
- Rate simulation 
Lietal. 201 lb E 
Li et al. 201 lb SO 



1<T 10" 

Total galaxy stellar mass pVL,] 

Figure 10. SN la rates per unit stellar mass, as a function of 
total galaxy stellar mass. The mass-normalized SN la rates for 
all galaxies in each mass bin are shown as black squares, while 
the rates for star-forming (passive) galaxies are shown as blue cir- 
cles (red diamonds). Vertical error bars are based on the Poisson 
uncertainty on the number of SNe la in the specific mass bin, 
and the horizontal error bars denote the range within which 68 
per cent of the galaxies fall within the mass bin. The solid curve 
shows the best-fitting SN la rate as a function of s tellar mass, as 
derived from a combination of a t^ 1 DTD and the lGallazzi et aD 
(2005) relation between mass and galaxy age. The shaded area 
is the confidence region resulting from the 68 per cent statistical 
uncertainty of the DTD amplitude The dashed (dot-dashed) 
curves are the Lll power-law fits to their mass-normalized SN la 
rates in local elliptical (SO) galaxies. 



of galaxy mass, for all galaxies in our sample, and separately 
for passive and star-forming galaxies. The mass-normalized 
SN la rates decrease with increasing galaxy mass. We thus 
confirm, with the present sample at z ~ 0.1, the similar 
dependence found by Lll for SNe la discovered in the local 
LOSS sample. 

|Kistler et al.| ( |2011| | have recently pointed out that the 
rate-mass relation found by Lll can be explained by the fact 
that massive galaxies formed earlier in the cosmic history, 



combined with a 



t~ 



DTD. We repeat here the Kistler 



et al. (20111 calculation, with some small differences. The 
mass-normalized SN la rate of galaxy i at cosmic time t 
can be expressed as the convolution of the SFH, S(t), and a 
DTD, ty(t), divided by the total stellar mass of the galaxy, 



M* , 



S(t')y(t-t')dt'. 



(6) 



Following Kistler et al. (20111, we use the Gallazzi et al 



(2005) relation between the stellar mass of a galaxy and its 
age to derive ages for a random selection of 10,000 galaxies 
with redshifts drawn from the redshift distribution of the 
SDSS DR7 sample. For each galaxy, we draw a galaxy age 
from a Gaussian distribution centred on the median values 



in table 2 from Gallazzi et al. ( 2005 1, with the 16/84 per cent 
values acting as the distribution's lower and upper standard 
deviations. Following Gallazz i et al.| ( [2005[ ), we use an expo- 
nential SFH of the form e~ a , with indices a drawn from a 
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Table 5. Mass-normalized SN la rates 



mass range 


median" 


SNc la 


SN la rate 6 


[10 10 M ] 


[10 10 M ] 




[SNuM] 


All galaxies 


< Af. < 3 
3 < Af. < 5.5 
5.5 < Af. < 10 
Af* > 10 

All masses 


1 3+ 1 ' 1 
J "°-1.Q 
4 1+0-9 

74 +i.7 

'•*-1.3 
20+ 24 


23 
22 
23 
22 
90 


171+U.U44 
ui ' i -0.035 
n 109+0. 035 
u.loz_ 02g 

096+ - 024 

n57+ - 015 

U.UB|-_ Q12 
0.10 ±0.01 




Star-forming g 


alaxies 




< M* < 2 
2 s£ Af, < 4 
4 s£ Af * < 8 

Af. 5= 8 
All masses 


U - 8 -0.6 

5-61!;? 

12±| 
,+6 

J -2 


11 
13 
13 
10 

47 


n 99 + O.O9 
u - z -0.06 

n 1S+006 
11+ 004 

u - il -0.03 
0fi1+ 026 

U.UD1_q Q19 

n 1 1Q+0-020 

U - il8 -0.017 


Passive galaxies 


< Af. < 4 
4 < Af » < 8 
8 < Af. < 12 
Af. > 12 
All masses 


9C+25 

111 2 . 2 


10 
10 
10 

13 
43 


u - lo -0.04 

076+ - 032 
u.u(o_ 024 

10+ 002 
u - lu -0.03 

0.0621^0?? 
0.082±g;g« 



a Uncertainties represent the mass region occupied 
by 68 per cent of the galaxies in that bin. 
b Rate uncertainties are Poisson uncertainties on 
the number of SNe la in each mass bin. 



flat distribution between and 1. The SFH is scaled to pro- 
duce the galaxy's formed mass, M/, over the period of time 
between a galaxy's formation time, t g , and t. We assume a 
power-law DTD with index —1, with the amplitude ^1 at 
t = 1 Gyr left as a free parameter. Equation[6]thus becomes: 



Mi 

M, 



*i 



'(t-t'^dt'. 



(7) 



As our measurements cover mainly old galaxies in the mass 
range ~ 10 9 - 10 12 M , we set Mf/M, = 2.3 ( |Bruzual fc 
|Charlot||2003[ ). The best-fitting DTD amplitude, with a re- 
duced is *i = 0.070 ± 0.016 x 10~ 12 Mq" 1 yr _ \ 
where the uncertainty of is the 68 per cent confidence 
region, defined as the range of $1 values within Ax 2 = ±1 
of the minimal \ 2 value. The form of the mass-normalized 
SN la rate as a function of mass is similar to the one ob- 
tained by |Kistler et al.| ( |2011[ ) , though with a shallower de- 
cline at the high-mass end. Fig. [To] shows the result of this 
simulation. 

Integrating the DTD with the best-fitting value for \ti 
over a Hubble time yields a number of SNe la per formed 
stellar mass of N/M» = 0.4 ± 0.1 x 10" 3 SNe M _1 . This 
value is lower than those obtained by Mil and M12, but 



consistent with the values obtained by Graur et al. (20111 
and P12. 



5.3 The Type la supernova volumetric rate 

We can convert the mass-normalized SN la rate, averaged 
over all masses and redshifts in our sample, to a volumetric 
rate, by multiplying it by the total cosmic mass density, as 
inferred from the galaxy stellar mass function (GSMF) mea- 
sured at z < 0.06 by Baldry et al. (20121. However, as the 



SDSS was a galaxy-targeted survey, we must account for its 
particular distribution of galaxy mass, M, shown in Fig. [I] 



weighted by the rate-mass relation we have found for these 
galaxies. We do this by multiplying our mass-normalized 
rate by the ratio of the integrated cosmic GSMF, B(M), 
to the SDSS DR7 galaxy-mass distribution, D(M) (which 
is normalized so that J D{M)MdM — 1), where both mass 
functions are weighted by the rate-mass relation R(M) de- 
scribed in Section \5. 21 



Ru 



Rla,M 



J B(M)R(M)MdM 
J D(M)R(M)MdM ' 



(8) 



The volumetric SN la rate, at a median redshift of 
z = 0.11, is J?i a ,v = 0.25±0.03 x 10~ 4 SNe yr _1 Mpc" 3 . The 
uncertainty in the rate derives from the Poisson uncertainty 
from the size of the SN la sample, along with a system- 
atic uncertainty of —1 SN la, as derived from the purity of 
our discovery and classification method, as detailed in Sec- 
tion [33] Varying the GSMF according to the uncertainties 



reported in table 1 of Baldry et al. (20121 adds a statistical 
uncertainty an order of magnitude smaller than the Poisson 
uncertainty, and is thus negligible. With an uncertainty of 
12 per cent, this is the most precise rate measurement at this 
redshift. The uncertainty, however, does not include the sys- 
tematic uncertainties in the mass estimates of the individual 
galaxies by T09, an uncertainty that is difficult to estimate. 
Our volumetric SN la rate measurement, along with mea- 
surements from other SN surveys in this redshift range, are 
shown in the upper panel of Fig. |11| Our volumetric rate 
is in excellent agreement with most previous measurements 
(e.g., M03; |Dilday et al.|2010[ ), but is more precise. Further- 
more, our rate merges smoothly with the extrapolation to 
low redshift of the rates measured by P12 (whose own low-z 
measurement has large errors). This is shown in the lower 
panel of Fig. |ll| Our low-2 rate measurement, together with 
the P12 measurements, thus provides now a precise picture 
of the SN la rate evolution from z = to 1. Following |Graur| 



et al. (20111 and P12, we convolve a power-law DTD with 
different cosmic SFHs (see Graur et al.|2011 section 8) and 
fit our volumetric SN la rate, together with those of P12 and 



Graur et al. (2011 see Fig. 111. We find a best-fitting in- 



dex of —1 ±0.1 ±0.16, where the uncertainties are statistical 
and systematic, respectively, and the systematic uncertainty 
derives from the different SFHs used in the fit. Taking into 
account the statistical and systematic uncertainties derived 



from the cosmic SFHs used in |Graur et al.| ( |201l[ ), N/M* 
lies in the range N/M* = (0.4-1.2) x 10 -3 SNe Mg" 1 . 

Our rate calculation is improved compared to the 
SN la rates derived by the two previous studies that have 
searched for SNe in SDSS spectra. M03 estimated the B- 
band- luminosity- weighted SN la rate at z w 0.1 using the 
average galaxy luminosity value of their sample, and a 
fixed visibility time estimated at t v — 20 days. They ob- 



tained 0.4 ± 0.2 h 2 SNuB (10 
100 km s" 1 Mpc -1 



si) y f 



where h 



Horesh et al. ( 2008 \ converted this rate 
into a volumetric rate of 0.24 ±0.12 x 10~ 4 SNe yr" 1 M pc~ 3 
using the redshift-dependent Botticcll a et al.| ( |2008[ ) lumi- 
nosity density function, but without taking into account the 
then-unknown rate-mass relation, which weights the par- 
ticular luminosity distribution of the galaxies in the SDSS 
DR1 sample. Similarly, Kll obtained their volumetric rate, 

0-569lS:2!s;±o.SS x 10 " 4 SNe y r_1 Mpc" 3 (the errors are 
statistical and systematic, respectively) by converting their 



B-band-luminosity- weighted rate of 0.472 



+0.081, +0.048 
-0.071,-0.039 



SNuB 



Supernovae in SDSS DR7 Spectra 15 



using the same redshift-dependent luminosity density func- 
tion, again without accounting for the effect of the rate-mass 
relation. Finally, neither M03 or Kf f made use of the de- 
tailed spectral information available for every galaxy in the 
sample, which can give the stellar masses of the galaxies 
(and which we have used for calculating rates, above) and 
the SFHs of the galaxies (which we use for recovering the 
DTD, below). 



5.4 The Type la supernova delay-time 
distribution 

The vbspa SFHs that are available for each of the individual 
galaxies that were effectively 'monitored' by our SN survey 
permit applying to the data the Mil algorithm for recover- 
ing the DTD. In this approach, briefly, Eq.[6]is approximated 
by a discretised form, 



N 



(9) 



where Ni is the expectation value of the SNe la in the ith 
galaxy, rriij is the stellar mass formed in the galaxy over the 
jth time interval, Vl" is a binned version of the DTD, and Ati 
is the visibility time of the ith galaxy. The sum is over the 
three vespa time bins: < 0.42 Gyr ('prompt'), 0.42-2.4 Gyr 
('intermediate'), and > 2.4 Gyr ('delayed'). The values of j 
are treated as free parameters that are scanned so as to best 
match Ni to the Poisson statistics of the observed sample, 
which has zero SNe in most of the monitored galaxies, and 
one SN in a small fraction of the galaxies. Mil, analysing 
the LOSS sample of Lll, detected significant prompt and 
delayed components to the DTD. Applying the method to 



a sample of SNe la from the SDSS-II survey ( Dilday et al. 
|20T0} [Sako et al.poTTj ), M12 obtained 4ct detections of all 
three DTD components. 

For the current SN la sample and the VESPA-derived 
SFHs for our galaxies, we obtain a strong detection of a 
delayed component in the age bin > 2.4 Gyr, of 4.5 ± 0.6 x 
10 -14 SNe Mq -1 yr _1 . Our recovered delayed component 
of the SN la DTD, along with measurements from other SN 
surveys, are shown in the Fig. |12| Our delayed component 
measurement is 2-3 times higher than the corresponding 
value in M12, but similar to the value obtained by Mil. It is 
possible that this result is affected systematically by the 'flux 
flow' problem we have identified in galaxies that hosted SNe. 
Some of the blue light of a host galaxy, light from a young 
or intermediate-age stellar population, is removed from the 
galaxy spectrum by the SN la template subtraction. As a 
result, some of the DTD signal could 'drift' from the prompt 
and intermediate bins to the delayed bin. 

We do not detect prompt and intermediate components 
to the DTD. In retrospect, this is unsurprising, since our 
galaxy sample is dominate by old and massive galaxies, with 
little ongoing star formation. 90 SNe la were hosted by the 
707,792 galaxies in our sample. Of these, only 5,394 galaxies 
have log(sSFR/yr _1 ) > —9.5, none of which hosted a SN la. 
The SDSS-II sample used in M12, on the other hand, is com- 
posed of 66,400 galaxies, which hosted a total of 132 SNe la. 
3,867 galaxies had log(sSFR/yr _1 ) > -9.5, seven of which 
hosted a SN. Comparing between the two samples, we would 
have expected (66,400/707,792) x (5,394/3,867) x (90/132) x 
7 = 0.6 SNe la in the highly-star-forming galaxies in our 
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Figure 11. SN la volumetric rate at z = 0.11 (red square) com- 
pared to rates from the literature. Upper panel: circles denote 
results by [C appellaro, Evans, & Turatto (1999); Hardin et al 



(2000) 



1 2010 i; 



Blanc et al.| (|2004|; |Horesh et al] ( |2008| l ; ]Dilday et al. 



Rodney & Tonry| | |26lO[ |; Lll; and P12. The M03 mea- 



surement, based on a spectroscopic SN search in SDSS DR1, 
the black diamond, and the Kll measurement, based on SDSS 
DR5, is the black upright triangle. The curve is the fit by P12 
to their measured SN la rates between z = 0.1 and 1. All ver- 
tical error bars include statistical and systematic uncertainties. 
The horizontal error bar indicates the range that encompasses 68 
per cent of the SDSS DR7 galaxies around the median redshift of 
z = 0.11. Lower panel: our low-z SN la rate (red square) compared 
to recent high-z r ates: |Dahlen, Strolger, fc Riess] | |2008[ l (down- 
turned triang les) , |Rodney fc Tonry||2 010) (upturned triangles), 
|Graur et al] ( |2011] (blue circles), |Barbary et al] ( |2012a| ) (left- 
facing triangles), and P12 (green diamonds). The solid curve is 
the same as above, and the dashed curve is the best-fitting SN la 
rate evolution derived from convolving the |Yiiksel et al] ( |2008[ ) 
cosmic SFH with a power-law DTD and fitting our volumetric 
rate measurement, together with those of P12 and |Graur et al] 
| |2011[ l . The most precise measurements in each redshift range are 
marked with filled symbols. 
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Figure 12. Recovered 'delayed' SN la DTD component (red 
square) compared to previous measurements from the literature: 
|Totani et al.| ( |2008| ) DTD measurements from SNe la discovered 
in 0.4 < z < 1.2 elliptical galaxies (diamonds); galaxy cluster 
surveys (triangles; |Maoz et al.|2010] and references therein; |Sand| 
|et al.|20 12); and the 'prompt', 'intermediate', and 'delayed' com- 
ponents recovered with the Mil method (purple; Mil and green; 
M12) circles. The solid curve is a i -1 DTD that, when convolved 
with the |Yiiksel et a.l.| ( |2008] ) cosmic SFH, provides the best fit to 
our volumetric rate, together with those of P12 and [Graur et al. 
120111 see Fig. [II}. The dashed curves are the systematic uncer- 
tainty on the amplitude of the DTD, derived from using different 
cosmic SFHs. 



sample, which is consistent with our having found none. 
With no detected SNe la in young galaxies, our method 
is unable to reconstruct the short- and intermediate-delay 
DTD components. 

For the time-integrated number of SNe la per formed 
stellar mass, the Mil recovery method gives N/M, — 
0.5 ±0.2 x 10~ 3 SNe M _1 . This is several times lower than 
the values found by Mil and M12, which is not surpris- 
ing, given that our survey misses the dominant prompt- and 
intermediate-delay SN la population. 



6 SUMMARY AND CONCLUSIONS 

We have surveyed the ~ 700,000 galaxy spectra in SDSS 
DR7 that have VESPA-derived SFHs, and discovered and 
classified 90 SNe la and 10 SNe IIP. The purity of our detec- 
tion and classification routine is such that there is, at most, a 
1 per cent systematic uncertainty in the classification of the 
SNe la. Using our SN la sample, together with the vbspa- 
derived stellar masses of the galaxies we survey, we measure 
the mass-normalized SN la rates at a median redshift of 
z = 0.11 and find that they follow the so-called 'rate-size', 
or rate-mass, relation, which was previously observed in the 
local Universe by LIT We confirm that more massive galax- 
ies have lower SN la rates per unit stellar mass, and that, for 
galaxies of a given total mass, the mass-normalized SN la 
rate is weakly dependent on the galaxy's sSFR. Following 
|Kistler et al.|p011[ ), we also confirm that this rate-mass re- 
lation is essentially a coincidental 'conspiracy'. On the one 



hand, less massive galaxies tend to form later in the history 
of the Universe ('downsizing'). On the other hand, numer- 
ous recent SN la rate studies point to a DTD that peaks at 
short delays and decreases monotonically with time ( |Maoz| 



|et al. | |2010[ |Maoz fc Badenes|[2fjT0l |Graur et al.|2011| Mil 



M12; |Barbary et al.|2012b||Badenes fc Maoz|2012| P12). The 

rate-mass relation results from the fact that, for more mas- 
sive galaxies, we are probing further down along the delayed 
tail of the DTD. Our mass-normalized SN la rates cover a 
mass range of ~ 10 9 -10 12 M . In the future, it would be 
interesting to extend the mass-normalized SN la rates to 
both lower and higher galaxy masses, where our (albeit sim- 
plified) calculation predicts that the rate per unit mass will 
plateau. 

We have also derived a mass-normalized SN la rate, 
averaged over all masses and redshifts in our galaxy sample, 
0.10 ±0.01 SNuM, and a volumetric rate, i?i a ,v(z = 0.11) = 
0.25 ± 0.03 x 10~ 4 SNe yr^ 1 Mpc" 3 . The latter is the most 
precise rate measurement at this redshift, and is consistent 
both with previous measurements from other surveys, and 
with the low-redshift extrapolation of the power-law fit to 
the SN la rate measurements of P12. 

As the SDSS DR7 galaxy sample is composed mainly 
of massive, old galaxies, applying the direct DTD recov- 
ery method of Mil yields only a 'delayed' component of 
4.5 ± 0.6 x 10~ 14 SNe M Q _1 yr" 1 in the delay-time range 
> 2.4 Gyr. This value is 2-3 times higher than the corre- 
sponding value in M12, but similar to the value obtained 
by Mil. The differences may be due to a systematic flux 
flow problem that exists in our SN discovery process. We 
find that the time-integrated number of SNe la per unit 
formed stellar mass derived from a power-law DTD with 
index —1 and scaled to yield the above value, is N/M* = 
0.5 ± 0.2 x 10~ 3 SNe M B "' , consistent with the values ob- 
tained by Graur et al. (2011 1 and P12, but lower than those 



obtained by Mil and M12. However, the value we find here 
is likely underestimated due to the fact that our survey mon- 
itors, as noted, a predominantly old and massive galaxy pop- 
ulation, and therefore samples only a fraction of the SNe la 
in the DTD. 

In terms of the progenitor question, our results lend fur- 
ther support to previous indications for a power-law DTD 
with an index ~ — 1. This support comes through the rate- 
mass relation that we measure (which is naturally explained 
with such a DTD), through the delayed DTD component we 
recover directly (this component is part of the above DTD), 
and through the precise volumetric rate we find (which 
merges smoothly with rates at higher redshifts and, com- 
pared to the cosmic SFH, implies such a DTD). It remains 
to be seen whether theoretical SD models may produce such 
a DTD, or whether this result indicates that many, most, or 
all SNe la arise through the DD channel. Further studies, 
using both SN la rates and other relevant techniques, will 
shed light on this issue. 

While large photometric SN surveys, such as the SDSS 
SN Survey ( |Dilday et al.|2010| [Sako et al.|2011[ ), Supernova 
Legacy Survey JNeill et al. 2006 P12), and the ongoing Palo- 
mar Transient Factory ( |Rau et al.|2009[ ) have generated large 
samples of SNe of all types, they require costly spectroscopic 
followup in order to robustly classify their candidates. Mas- 
sive spectroscopic surveys, such as previous data releases of 
the SDSS or the ongoing SDSS-III ( |Eisenstein etllpnTj ) , 
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provide an ideal platform for our code, with which one can 
assemble large transient samples at no extra cost. The de- 
tection and classification method introduced here can be 
improved upon by, e.g., optimizing it for other types of SNe 
and by replacing the template library with bases of SN eigen- 
spectra, similar to the salt2 templates. Furthermore, we 
will search for ways to discover new types of transients that 
cannot be fitted with current spectral templates. Finally, al- 
though the current work has been an archival search, there is 
no obstacle to discovering SNe in spectra in real time, while 
the SNe can still be followed up. We are currently apply- 
ing such an approach to spectra from the SDSS-III Baryon 
Oscillation Spectroscopic Survey ( Dawson et al.|2012 l. 
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